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Laser  Induced  Fluorescence  Measurements  within  a 
Laboratory  Hall  Thruster 

W.A.  Hargus,  Jr.*  and  M.A.  Cappelli 
Mechanical  Engineering  Department 
Thermosciences  Division 
Stanford  University 
Stanford,  CA  94305 


Abstract 

In  this  paper,  we  describe  the  results  of  a  study  of  laser  induced  fluorescence  velocimetry  of  ionic 
xenon  in  the  plume  and  interior  acceleration  channel  of  a  laboratory  Hall  type  thruster  operating  at  powers 
ranging  from  250  to  725  W.  Optical  access  to  the  interior  of  the  Hall  thruster  is  provided  by  a  1  mm  axial  slot 
in  the  insulator  outer  wall.  Axial  ion  velocity  profiles  for  four  discharge  voltages  (100  V,  160  V,  200  V,  250 
V)  are  measured  as  are  radial  velocity  profiles  in  the  near  field  plume.  Internal  neutral  xenon  axial  velocity 
profiles  are  also  measured  at  these  conditions.  For  each  test  condition,  the  implied  axial  electric  field  is  cal¬ 
culated  from  the  measured  ionic  velocity  profiles.  These  results  are  compared  to  previous  plasma  potential 
measurements  performed  with  an  emissive  probe.  The  correspondence  between  the  two  sets  of  measure¬ 
ments  indicates  that  the  ionization  and  acceleration  regions  within  the  thruster  are  separated  to  some  degree. 


Introduction 

Due  to  their  high  specific  impulse  and  high 
thrust  efficiencies.  Hall  thrusters  are  now  being 
considered  for  use  on  commercial,  research,  and 
military  spacecraft.  This  technology  provides  eco¬ 
nomic  advantages  for  a  number  of  missions  and  its 
use  can  be  translated  into  lower  launch  mass,  longer 
time  on  station,  or  larger  payloads  [1 [. 

There  is  a  need  for  increased  understand¬ 
ing  of  the  complex  phenomena  that  govern  the 
operation  of  Hall  thrusters.  In  order  to  more  fully 
understand  the  physics  in  these  discharges,  several 
laboratory  model  Hall  thrusters  have  been  con¬ 
structed  at  Stanford  University.  These  thrusters 
have  served  as  test  articles  for  model  development 
and  advanced  plasma  diagnostics  including  laser 
induced  fluorescence  (LIF),  probes  of  various 
types,  and  thrust  measurements  [2-8], 

Laser  based  techniques  have  been  devel¬ 
oped  to  nonintrusively  probe  neutral  and  ionized 
xenon  [6-8].  Such  measurements  in  the  plumes  of 
Hall  and  other  types  of  ion  thrusters  provide  infor¬ 
mation  on  the  plasma  discharge  that  is  useful  in 
developing  insight  into  the  physical  processes 
occurring  within  these  devices.  Similar  optical 
diagnostic  measurements  have  been  previously 
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employed  to  examine  plasma  properties  in  other 
electric  propulsion  devices.  For  example,  the  hydro¬ 
gen  arcjet  has  been  extensively  studied  using  lasers 
to  measure  velocity,  temperature,  and  electron  num¬ 
ber  density  [9[.  The  high  spatial  resolution  of  single 
point  laser-induced  fluorescence  is  essential  in 
probing  nonuniform  plasma  environments  such  as 
those  in  Hall  thrusters  and  other  electric  propulsion 
devices. 

Theory 

Laser  Induced  Fluorescence 

The  interaction  of  a  laser  beam  with  a 
plasma  may  involve  the  optical  excitation  of  a  num¬ 
ber  of  atoms  to  a  higher  energy  state.  The  excitation 
is  more  likely  to  occur  if  the  laser  is  tuned  to  the 
energy  difference  hv12  between  an  upper  and  lower 
excitation  level.  The  interaction  can  be  investigated 
by  either  monitoring  the  resulting  reduction  in  laser 
power  following  propagation  through  the  plasma 
(absorption  process),  or  by  monitoring  the  subse¬ 
quent  spontaneous  emission  as  the  resulting  excited 
state  relaxes  to  a  lower  state  (laser-induced  fluores¬ 
cence,  LIF).  Monitoring  the  fluorescence  as  the 
laser  is  tuned  over  the  transition  provides  a  measure 
of  the  fluorescence  excitation  line  shape  and  is 
favored  for  the  higher  spatial  resolution  that  it 
affords  in  the  determination  of  plasma  parameters. 
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The  spatial  resolution  for  LIF  is  determined  by  the 
intersection  of  the  probe  laser  beam  with  the  optical 
collection  volume. 

If  an  absorber  has  a  velocity  component  u 
along  the  axis  of  the  laser  beam,  it  will  absorb  the 
light  at  a  frequency  shifted  from  that  of  stationary 
absorbers  due  to  the  Doppler  effect.  The  magnitude 
of  this  frequency  shift  depends  on  the  velocity 
along  the  laser  beam  axis  by. 


where  c  is  the  speed  of  light.  The  Doppler  shift  of  a 
species’  fluorescence  profile  away  from  the  line 
center  v12  of  stationary  absorbers  is  in  proportion 
to  the  species  bulk  velocity. 

The  measured  fluorescence  signal  is  given 
by  fl  0-1 2  J 

Sf  =  rld«chvl2A2t  N2  (2) 

where  q12  is  the  efficiency  of  the  detection  system, 
ac  accounts  for  factors  involving  the  collection 
system,  and  A21  is  the  Einstein  coefficient  for 
spontaneous  emission  of  the  relevant  transition.  For 
low  laser  intensities,  rate  equation  analysis  indi¬ 
cates  that  the  upper  level  population  N2 ,  and  there¬ 
fore  the  fluorescence  signal,  is  linearly  dependent 
on  laser  intensity  at  steady  state,  i.e., 

N2  -  IVB12<(>V  (3) 

where  lv  is  the  spectral  irradiance  at  frequency  v , 
B]2  is  the  Einstein  stimulated  absorption  coeffi¬ 
cient,  and  4>v  is  the  transition’s  spectral  line  shape 
which  accounts  for  the  variation  of  the  absorption 
or  laser  excitation  with  frequency.  The  line  shape  is 
determined  by  the  environment  of  the  absorbing 
atoms,  so  an  accurate  measurement  of  the  line 
shape  function  can  lead  to  the  determination  of  var¬ 
ious  plasma  parameters.  However,  for  velocity  mea¬ 
surements,  partially  saturated  fluorescence  with  a 
distorted  line  shape  is  still  capable  of  providing  a 
reasonable  measure  of  the  mean  ion  velocity.  This 
was  experimentally  verified  with  two  saturation 
studies  which  also  determined  velocity.  The  varia¬ 
tion  of  the  velocities  were  found  to  be  less  than  the 
experimental  uncertainty  for  the  ions  (±500  m/s),  or 
for  the  neutrals  (±60  m/s). 


Several  factors  can  affect  the  line  shape 
and  give  rise  to  a  broadening  and/or  shift  of  the 
spectral  line.  In  high  temperature  plasmas,  the  most 
significant  of  these  is  due  to  the  Doppler  effect  and 
gives  rise  to  Doppler  broadening.  Doppler  broaden¬ 
ing  reflects  the  fact  that  the  absorbing  species  is  dis¬ 
tributed  in  velocity  according  to  the  kinetic 
temperature,  Tkin .  When  the  absorbing  species  is 
characterized  by  a  Maxwellian  velocity  distribu¬ 
tion,  Doppler  broadening  results  in  a  Gaussian  line 
shape.  Collisional  interactions  between  the  absorb¬ 
ers  and  other  species  in  the  plasma  give  rise  to  spec¬ 
tral  line  shapes  that  are  often  Lorentzian.  This 
includes  interactions  with  charged  particles  (Stark 
broadening)  and  uncharged  particles  (van  der  Waals 
broadening).  If  both  Doppler  broadening  and  colli¬ 
sion  broadening  are  important  and  independent,  the 
resulting  line  shape  is  a  convolution  of  the  Gaussian 
and  Lorentzian  line  shape  into  a  Voigt  line  shape 
U2J. 

The  line  shape  is  an  intrinsic  property  of 
the  absorbers,  whereas  the  fluorescence  excitation 
line  shape  is  the  variation  in  the  detected  fluores¬ 
cence  signal  with  frequency  as  the  laser  is  tuned 
across  the  absorption  line  feature.  If  the  laser  exci¬ 
tation  significantly  perturbs  the  populations  in  the 
coupled  levels,  it  is  said  to  be  saturating  the  transi¬ 
tion  and  the  fluorescence  signal  is  then  a  nonlinear 
function  of  laser  intensity.  In  cases  where  the  laser 
intensity  is  significantly  below  the  saturation  level 
and  the  laser  linewidth  is  small  compared  to  the 
measured  linewidth,  the  fluorescence  excitation  line 
shape  reflects  the  spectral  absorption  line  shape  as 
given  by  Eqns.  2  and  3.  When  the  laser  intensity  is 
sufficiently  high  to  saturate  the  transition,  the  fluo¬ 
rescence  excitation  line  shape  is  broader  than  the 
spectral  line  shape  and  the  fluorescence  intensity  is 
less  than  it  would  be  if  it  were  linear  with  the  laser 
intensity  Iv .  The  saturation  intensity,  defined  as  that 
intensity  which  produces  a  fluorescence  signal  half 
of  what  it  would  be  if  the  fluorescence  was  linear 
with  Iv,  depends  inversely  on  the  line  strength  of 
the  particular  line.  Stronger  transitions  have  a 
smaller  saturation  intensity  and  thus  a  larger  satura¬ 
tion  effect  for  a  given  laser  intensity. 

Hvperfine  Structure 

The  seven  xenon  isotopes  each  have  a 
slight  difference  in  their  transition  energies  due  to 
their  differences  in  mass.  The  odd  mass  isotopes  are 
further  spin  split  due  to  nuclear  magnetic  dipole  and 
nuclear  electric  quadrapole  moments  [131. 
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Many  nuclei  which  have  an  odd  number  of 
protons  and/or  an  odd  number  of  neutrons  posses  an 
intrinsic  nuclear  spin  Ih/2 n .  Where  I  is  integral, 
or  half-integral,  dependent  on  if  the  atomic  mass  is 
even  or  odd,  respectively  [14],  For  nuclei  with 
nuclear  spin,  there  exists  an  interaction  of  the 
nucleus  with  the  electron  shell.  The  interaction 
leads  to  the  splitting  of  levels  with  angular  momen¬ 
tum  J  into  a  number  of  components,  each  corre¬ 
sponding  a  specific  value  of  the  total  angular 
momentum  F  [15]. 

F  =  I  +  J  (4) 

As  a  result  of  this  interaction,  F  is  a  conserved 
quantity  while  I  and  J  are  not.  The  interaction  is 
very  small  allowing  the  hyperfine  splitting  of  each 
level  to  be  taken  independently  of  the  other  levels. 
The  number  of  nuclear  spin  split  hyperfine  compo¬ 
nents  is  21  +  1  if  J  2 1  and  2J  +  1  if  J  <  I  with  F  val¬ 
ues  shown  below  [15,16], 

F  =  J  +  l,  J  +  I-l, ...,  |J-I|  (5) 

In  turn,  the  selection  rules  for  F  are  given 

below. 

AF  =  0,  ±  1  (6) 

AF  *  0 ,  if  F  =  0  (7) 

Two  constants  are  associated  with  hyper¬ 
fine  nuclear  spin  splitting.  These  are  the  A  hyper¬ 
fine  structure  constant  which  represents  the  nuclear 
magnetic  dipole  effect  on  the  atom,  and  the  B 
hyperfine  structure  constant  which  is  associated 
with  the  nuclear  electric  quadrapoie  moment  of  the 
atom  which  will  only  be  present  if  I  2  1  .  The  rela¬ 
tive  energy  of  the  spin  split  states  depends  on  the 
sign  of  A .  In  atoms  with  A  >  0 ,  the  state  with  the 
highest  value  of  F  has  the  highest  energy,  while  for 
atoms  with  A  <  0 ,  the  state  with  the  lowest  value  of 
F  has  the  highest  energy  1 16].  The  energy  level 
shift  AEm  associated  with  the  magnetic  dipole  of 
the  nucleus  is  given  by  1 14] 

AEm(F)  =  IA[F(F+1)-J(J+1)-I(I+1)]= 

(8) 

Additionally,  the  energy  spacing  between  succes¬ 
sive  levels  F-  1  and  F  may  be  shown  to  be  propor¬ 
tional  to  F . 


AEm(F)-AEm(F-1)  =  AF  (9) 

If  I  2  1  ,  the  nucleus  will  have  an  electric  quadra¬ 
poie  moment  and  a  related  hyperfine  splitting  con¬ 
stant  B  which  produces  an  additional  hyperfine 
splitting  with  energy  linear  in  C(C+  1)  where  C  is 
defined  in  Eqn.  8. 

AEF  =  AEm  +  AEq  =  (10) 

+  b[c(C  +  l)-|j(J+  1)1(1+  1)] 

Where  AEF  is  the  combined  nuclear  spin  split 
energy  level  shift  combining  the  effect  from  the 
nuclear  magnetic  dipole  moment,  AEM ,  and  the 
effect  of  the  electric  quadrapoie  moment,  AEQ  [  14]. 

It  should  be  noted  that  the  center  of  gravity  of  the 
hyperfine  levels  lies  at  a  position  of  the  unsplit  level 
j  [151. 

(2F  +  1  )AEf  =  0  (11) 

Due  to  close  energy  spacing  of  nuclear 
spin  split  levels,  near  ideal  coupling  between  I  and 
J  occurs  in  most  hyperfine  structure.  Therefore,  the 
relative  intensity  rules  derived  in  the  literature  for 
Russel-Saunders  coupling  are  appropriate  for 
hyperfine  splitting  [  16]. 

For  where  J  -  1  -»  J , 

F-  1  -»  F : 

(J  +  F+I+  1)(J  +  F+I)(J  +  F-I)(J  +  F-I-1) 

F 

F-*F: 

(J  +  F  +  I+  1)(J  +  F-I)(J-F  +  1)(J-F-1-I)(2F+  1) 
F(F+  1 ) 

F+  I  -  F: 

(J  —  F+1)(J  —  F+l  —  1)(J  —  F  —  I  —  1)(J  —  F  —  I  —  2) 

(F  +  1) 

For  where  J  — *■  J , 

F  -  1  -*  F : 

(J  +  F+1+  !)(JtF-l)(J-FtI  +  I )( J  —  F  —  I ) 

F 

F  — F: 

[J(J+  1 )  +  F(F  +  1)  +  I(1+  1 )  ]2(  2F  +  1) 

F(F+  1) 
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( J  +  F  +  I  -  2)(  J  +  F - 1  +  1 )( J  -  F  +  I)( J  -  F  -  I  -  1 ) 
(F+l) 

(12) 

The  relative  intensities  of  the  isotope 
shifted  transitions  are  proportional  to  each  isotope’s 
relative  abundance  shown  in  Table  1.  The  relative 
intensities  of  the  nuclear  spin  split  hyperfine  split¬ 
ting  are  governed  by  two  summation  rules  [15]. 
First,  the  sum  of  the  intensities  all  the  lines  of  the 
hyperfine  structure  of  a  transition  J  -*  J’  originating 
from  a  component  F  of  the  level  J  is  proportional 
to  the  statistical  weight  of  this  component,  2F  +  1 . 
Second,  the  sum  of  the  intensities  of  all  the  lines  of 
the  hyperfine  structure  the  transition  J  -*  J'  ending 
on  the  component  F  of  the  level  J’  is  proportional 
to  the  statistical  weight  of  this  component, 
(2F  +  1 ) .  With  these  two  sum  rules,  a  system  of  lin¬ 
ear  equations  are  solved  for  the  relative  intensities 
of  the  nuclear  spin  split  components  of  each  iso¬ 
tope. 

Xenon  Spectroscopy 

The  isotopic  and  nuclear-spin  effects  con¬ 
tributing  to  the  hyperfine  structure  of  the 
6s[3/2]^  -6p[3/2]2  neutral  xenon  transition  at 
823.2  nm  produce  21  individual  lines.  Similarly,  the 
5d[4]7/2-6p[3]5/2  ionic  xenon  transition  at  834.7 
nm  shown  in  Fig.  1  has  a  total  of  19  isotopic  and 
spin  split  components.  The  hyperfine  splitting  con¬ 
stants  are  only  known  for  a  limited  set  of  levels  sev¬ 
eral  of  which  are  shown  in  Table  2.  The  neutral 
823.2  nm  transition  has  all  of  these  measured  and 
tabulated  in  the  literature.  The  ionic  834.7  nm  tran¬ 
sition  only  has  data  on  the  nuclear  spin  splitting  of 
the  upper  state  and  no  available  information  on  the 
transition  dependent  isotope  shifts. 

For  LIF  measurements  primarily  aimed  at 
determining  velocities  within  the  plasma  flow,  it  is 
often  convenient  to  probe  more  accessible  transi¬ 
tions  for  which  there  is  incomplete  knowledge  of 
the  isotopic  and  nuclear  spin  splitting  constants. 


Fig.  1.  Nuclear  spin  splitting  of  the 
5d[4]7/2  -  6p[3]5/2  transition  of  singly  ionized 
xenon  at  834.7  nm. 


Manzella  has  shown  that  the  ionic  xenon  transition 
at  834.7  nm  can  be  used  to  make  velocity  measure¬ 
ments  in  the  plume  of  a  Hall  thruster  [17,18],  An 
additional  convenient  feature  of  this  transition  is  a 
strong  line  emanating  from  the  same  upper  state. 
This  transition  at  541.9  nm  allows  for  nonresonant 
fluorescence  collection. 


Table  1 :  Xenon  natural  isotope  abundances,  nuclear  spins, 
and  representative  isotope  shifts  for  Xe  I  and  Xe  II 
transitions  relative  to  132Xe  in  MHz  [  14,19-20 1. 


%  Natural 
Abundance 

Nuclear 

Spin 

Xe  I 

823.2nm 

Isotope 

Shift 

Xe  II 

605.  Inm 
Isotope 
Shift 

124Xe 

0.10 

0 

270  ±  12 

336.6  ±3.9 

126Xe 

0.09 

0 

189  ±  12 

252.4  ±3.9 

,28Xe 

1.91 

0 

1 17  ±  12 

172.0  ±6.8 

129Xe 

26.4 

1/2 

123±  18 

113.7  ±5.0 

,30Xe 

4.1 

0 

54  ±  12 

83.6  ±4.5 

13lXe 

21.2 

3/2 

81  ±  18 

16.7  ±5.0 

l32Xe 

26.9 

0 

0 

0 

I34Xe 

10.4 

0 

-30  ±12 

-75.8  ±4.0 

l36Xe 

8.9 

0 

-147  ± 

12 

-140.9  ± 

7.3 

Table  2:  Hyperfine  spin  splitting  constants  for  Xe  I  and  Xe 
II  in  MHz  123-25). 


State 

Al» 

AI3I 

B131 

Xe  1 

6s[3/2]£ 

-2384  ±2 

707  ±  1 

525  ±2 

Xe  I 

6p[3/2]2 

887  ±  1 

265  ±3 

30  ±3 

Xe  11 

5d[3]7/2 

-502  ±4 

150  ±  2 

70  ±  13 

Xe  11 

6p[2]5/2 

-1633  ±  It 

487  ±  4 

-129  ±  16 

Xe  11 

6p[3]7/2 

-1387  ±9 

409  ±2 

-117±  10 

Xe  11 

6p[3]5/2 

-1 164  ±  5 

346  ±2 

-230  ±  1 1 
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Test  Apparatus 

Hall  Thruster 

The  thruster  used  in  this  study  has  been 
described  extensively  elsewhere  [7].  It  has  4  outer 
magnetic  windings  consisting  of  89  mm  long,  25 
mm  diameter  cores  of  commercially  pure  iron 
wrapped  with  6  layers  of  22  gauge  insulated  copper 
magnet  wire.  The  inner  core  is  also  25  mm  in  diam¬ 
eter  and  89  mm  in  height  and  has  12  layers  of  mag¬ 
net  wire.  The  depth  of  the  electrical  insulator  is  84 
mm.  The  insulator  is  constructed  from  2  sections  of 
cast  99.9%  alumina  tubing  cut  to  length.  These  two 
pieces  are  cemented  to  a  machinable  alumina  plate 
attached  to  the  back  plate  of  the  thruster  with  non¬ 
conducting  fasteners.  A  photograph  of  the  assem¬ 
bled  thruster  is  shown  in  Fig.  2. 

A  second  insulator  was  constructed  for  the 
modified  Flail  thruster.  It  is  identical  to  the  insulator 
described  above  with  a  slot  approximately  1  mm 
wide  along  the  length  of  the  outer  wall  of  cast  alu¬ 
mina  cut  with  a  diamond  saw  prior  to  assembly.  The 
slot  is  used  to  provide  optical  access  to  the  interior 
of  the  Hall  thruster.  Operation  of  the  Hall  thruster 
with  the  slotted  insulator  does  not  appear  to  differ 
significantly  from  the  operation  with  the  unslotted 
insulator.  The  current  voltage  characteristics  are 
identical  and  there  are  no  apparent  differences  in 
thruster  operation.  A  photograph  of  the  slotted  insu¬ 
lator  is  also  shown  in  Fig.  2.  Here,  the  insulator  is 
shown  after  testing  has  been  completed.  The  dark 
deposits  on  the  insulator  arc  the  result  of  the  partial 
carbonization  of  the  phenolic  thruster  mount. 


The  optical  access  provided  by  the  slot  is 
blocked  near  the  exit  plane  by  the  front  plate  of  the 
magnetic  circuit.  The  front  plate  is  not  cut  since  this 


Fig.  2.  Photograph  of  Hall  thruster  with  slotted 
insulator  beside  it. 


would  modify  the  local  magnetic  field.  The  mag¬ 
netic  field  is  one  of  the  most  important  parameters 
determining  the  behavior  of  the  plasma  in  the  Hall 
thruster  and  any  change  in  the  field  would  immedi¬ 
ately  change  the  characteristics  of  the  plasma  It  is 
felt  that  the  4  mm  field  of  view  blocked  by  the  front 
plate  of  the  magnetic  circuit  is  less  important  than 
preserving  the  uniformity  of  the  magnetic  field 
within  the  Hall  thruster. 

The  thruster  is  mounted  on  a  two  axis 
translation  system.  In  the  vertical,  the  thruster  has  a 
range  of  travel  of  approximately  30  cm.  In  the  axial, 
or  horizontal,  direction,  the  thruster  is  constrained 
to  a  total  travel  of  approximately  6  cm.  Both  stages 
have  resolutions  on  the  order  of  10  pm,  although 
the  repeatability  is  considerably  courser.  For  inter¬ 
nal  LIF  measurements  with  the  slotted  insulator,  the 
thruster  was  mounted  on  a  platform  supported  by 
the  axial  translation  stage  and  a  linear  ball  bearing 
pillow  block  which  provides  unimpeded  optical 
access. 

Laser  Induced  Fluorescence 

The  experimental  apparatus  used  for  the 
laser  induced  fluorescence  measurements  consists 
of  a  tunable  Coherent  899-21  single  frequency  tita¬ 
nium  sapphire  laser.  The  laser  is  actively  stabilized 
to  provide  line  widths  on  the  order  of  1  MHz  with 
near  zero  frequency  drift.  Scan  widths  of  up  to  20 
GHz  can  be  realized  at  center  wavelengths  between 
680  to  1060  nm.  The  titanium  sapphire  laser  is 
pumped  by  Coherent  solid  state  Verdi  pump  laser. 
The  pump  laser  provides  5  W  of  single  mode  pump 
power  at  532  nm.  The  laser  wavelength  was  moni¬ 
tored  by  a  Burliegh  Instruments  WA-1000  scanning 
Michaelson  interferometer  wavemeter  with  a  reso¬ 
lution  of  0.01  cm'1.  Figure  3  shows  a  schematic  of 
the  laser  configuration. 

The  probe  beam  is  directed  into  the  Hall 
thruster  plume  by  a  series  of  mirrors.  The  slightly 
divergent  beam  (1.7  milliradians  full  angle)  is 
focused  to  a  submillimeter  beam  waist  by  a  50  mm 
diameter,  1 .5  m  focal  length  lens.  For  radial  veloc¬ 
ity  measurements,  the  probe  beam  enters  through  a 
side  window  and  is  focused  by  a  50  mm  diameter, 
50  cm  focal  length  lens.  The  collection  optics  for 
both  radial  and  axial  velocity  measurements  consist 
of  a  75  mm  diameter,  60  cm  focal  length,  collimat¬ 
ing  lens.  The  collected  light  is  then  focused  on  to 
the  entrance  slit  of  a  0.5  m  Ebert-Fastie  mono¬ 
chrometer  with  a  50  mm  diameter,  30  cm  focal 
length,  lens.  An  optical  field  stop  is  placed  between 
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Rjr>g  *e  Gta*  Photodiode 


Fig.  3.  LIF  velocimetry  apparatus  with  laser  system 
and  collection  optics. 


Pump 

Beam 


Hall  Thruster  Translation  Stages 


PMT-,-  -fc-  Rxxsmg  Loos 

J?,t»;MonodTOn«tBr 
GaSi5SBER^SSSE8a,gaUit2 
Axial  LIF  Veioctmetry  Measurements 


Fig.  4.  Stanford  high  vacuum  facility  and  LIF  col¬ 
lection  optics  with  cutaway  view  of  Flail  thruster. 


the  two  lenses  to  match  the  F Iti  of  the  optical  train 
with  that  of  the  monochrometer.  A  schematic  of  the 
collection  optical  train  is  shown  in  Fig.  4.  The 
monochrometer  is  used  as  a  narrow  band  optical  fil¬ 
ter  so  that  only  light  from  the  transition  of  interest  is 
collected.  With  entrance  and  exit  slits  full  open  (425 
/mi),  the  600  groove/mm  plane  grating  blazed  for 
600  nm  within  the  monochrometer  allows  the 
Hamamatsu  R928  photomultiplier  tube  (PMT)  to 
sample  a  wavelength  interval  of  approximately  1 
nm.  The  orientation  of  the  monochrometer  allows 
the  height  of  the  slits  to  define  the  length  of  the 
probe  beam  along  which  the  fluorescence  is  col¬ 
lected.  The  sample  volume  for  all  data  presented  in 
this  work  is  approximately  100  /im  in  diameter  and 
2  mm  in  length.  Prior  to  each  test  and  with  the 
PMT  housing  removed,  a  HeNe  laser  is  reversed 
through  the  entire  collection  optical  train  to  insure 
the  alignment  with  the  probe  beam.  For  neutral  LIF 
velocimetry  measurements,  a  portion  of  the  probe 
beam  is  split  from  the  main  beam,  passed  through  a 
xenon  glow  discharge  tube,  and  used  as  an  station¬ 
ary  absorption  reference.  A  silicon  photodiode 
monitors  the  absorption  signal.  This  use  of  the  glow 
discharge  tube  is  only  possible  for  neutral  xenon. 
The  glow  discharge  does  not  support  a  sufficient 
population  of  excited  state  ions. 


The  LIF  signal  is  collected  using  a  Stan¬ 
ford  Research  Systems  SRS-850  digital  lockin 
amplifier.  The  probe  bean  is  chopped  by  an  SRS- 


540  optical  chopper.  The  absorption  signal  from  the 
stationary  reference  is  collected  using  an  SRS-530 
lockin  amplifier.  Data  from  the  absorption  signal, 
laser  power  output,  and  the  wavemeter  are  stored  on 
the  SRS-850  using  3  available  analog  inputs  along 
with  the  LIF  signal.  Typical  tests  consist  of  a  12-20 
GHz  scan  of  the  probe  laser  frequency  over  a  3 
minute  period.  The  beam  is  chopped  at  a  frequency 
of  1.5  kHz.  Both  lockin  amplifiers  use  1  s  time  con¬ 
stants.  Data  is  sampled  at  8  Hz,  producing  four 
traces  of  approximately  2,000  points  for  each  veloc¬ 
ity  data  point.  Several  unsaturated  traces  using 
lower  laser  intensities,  10  s  time  constants,  and  10 
min  scans  were  also  performed. 

Results  and  Analysis 

Operating  Conditions 

The  Hall  thruster  was  operated  at  four  con¬ 
ditions.  At  each  of  these  conditions,  the  peak  mag¬ 
netic  field  was  measured  to  be  approximately  125 
G,  the  mass  flow  to  the  anode  was  2  mg/s,  mass 
flow  to  the  cathode  was  0.3  mg/s.  The  test  condi¬ 
tions  corresponded  to  discharge  voltages  of  100, 

1 60,  200.  and  250  V.  The  anode  currents  for  these 
conditions  were  2.1,  2.4,  2.6,  and  2.9  A,  respec¬ 
tively.  The  total  power  consumed  by  the  cathode 
and  magnet  circuit  was  approximately  30  W.  It 
should  be  noted  that  the  power  dissipated  in  the  bal¬ 
last  resistors  on  the  anode  and  cathode  keeper  lines 
(~10  W)  are  not  included  in  these  calculations. 

Position  Reference 

All  spatially  resolved  measurements  are 
referenced  to  a  two  coordinate  system.  The  position 
in  the  radial  coordinate  is  referenced  to  the  differ¬ 
ence  of  the  radial  location  and  the  radial  location 
corresponding  to  the  acceleration  channel  center 
using  the  variable  D  which  is  defined  as  positive 
toward  the  thruster  centerline.  The  axial  coordinate 
is  given  by  Z  which  is  the  distance  from  the 
thruster  exit  plane  and  is  defined  as  positive  along 
the  thrust  vector. 

Saturation  Study 

The  5d[4]7/2-6p[3]5/2  ionic  xenon  tran¬ 
sition  at  834.7  nm  was  probed  to  extract  local 
velocity  data  from  the  doppler  shift  of  the  measured 
fluorescence.  Figure  5  shows  the  saturation  curve  of 
the  transition  at  a  location  13  mm  from  the  exit 
plane  and  at  the  center  of  the  acceleration  channel. 
Here,  the  laser  beam  is  propagating  normal  to  the 
thrust  vector  and  is  measuring  the  radial  velocity 
component  of  the  xenon  ions.  A  typical  saturated 
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Fig.  5.  Saturation  curve  of  radial  ionic  measure¬ 
ments  at  Z  =  13  mm  and  D  =  0  mm. 


trace  used  to  determine  the  velocity  in  the  probed 
volume  is  compared  to  an  unsaturated  trace  from 
the  same  location  in  Fig.  6.  Both  traces  are  normal¬ 
ized  to  unity  peak  signal. 

Velocity  measurements  for  each  of  the  sat¬ 
urated  line  shapes  making  up  Fig.  5  are  shown  in 
Fig.  7.  The  mean  velocity  for  this  data  points  yields 
a  value  of  -62  m/s  with  a  standard  deviation  of  65 
m/s  and  a  range  of  127  m/s.  The  absolute  accuracy 
of  the  measurements  is  determined  by  the  uncer¬ 
tainty  of  the  wavelength  meter  and  is  taken  to  be 
±500  m/s.  The  most  important  conclusion  that  can 
be  drawn  from  this  data  is  that  the 
5d[4]7/2-6p[3]5/2  ionic  transition  of  xenon  pro¬ 
vides  useful  LIF  velocimetry  data  when  partially 
saturated  which  is  within  the  uncertainties  of  the 
measurement  of  the  Doppler  shift.  This  allows  the 
collection  of  saturated  fluorescence  signals,  maxi- 


Fig.  6.  Saturated  and  unsaturated  ionic  LIF 
traces 


Fig.  7.  Invariance  of  measured  radial  ionic  velocity 
in  saturation  study. 


mizing  the  signal  to  noise  ratio  and/or  allowing  for 
faster  scans  to  be  taken  of  the  line  while  still  being 
able  to  extract  valuable  velocity  data  from  the 
acquired  fluorescence  signal.  Similar  results  were 
obtained  in  a  neutral  saturation  study. 

Ionic  LIF  Velocimetry 

The  axial  LIF  velocimetry  data  consists  of 
two  sets  of  measurements  taken  with  the  slotted 
insulator.  The  first  set  consists  of  ionic  velocity 
measurements  taken  externally  extending  from  the 
exit  plane  to  approximately  Z  =  35  mm.  The  second 
set  of  data  points  consist  of  internal  axial  velocity 
data  from  the  exit  plane  to  approximately  Z  =  -65 
mm.  These  limits  are  imposed  by  the  limited  range 
of  the  translation  stage  providing  axial  motion  of 
the  Hall  thruster.  Between  each  data  set,  the  transla¬ 
tion  stages  had  to  be  repositioned  so  that  the  subse¬ 
quent  data  set  could  be  taken.  Data  points  were 
taken  every  2.5  mm  with  a  sample  probe  volume 
100  pm  in  diameter  and  2  mm  long.  Overlap  of  the 
two  data  traces  ensure  that  the  measurements  are 
continuous.  For  several  of  the  test  conditions  exam¬ 
ined,  profiles  of  the  axial  velocity  across  the  coordi¬ 
nate  D  are  also  examined. 

The  complete  axial  velocity  profiles  for 
the  four  cases  examined  are  shown  in  Fig.  8.  The 
error  bars  correspond  to  the  uncertainty  associated 
with  the  determination  of  the  magnitude  of  the 
Doppler  shift  (500  m/s).  The  axial  velocity  profiles 
exhibit  a  typical  behavior.  The  velocity  is  near  zero 
near  the  anode  (Z  =  -78  mm),  and  begins  to  rise 
near  Z  =  -10  mm.  The  ions  are  rapidly  accelerated 
in  the  region  of  the  exit  plane  and  reach  their  full 
velocity  in  the  neighborhood  of  Z  =  20  mm.  This 
latter  position  corresponds  to  the  location  of  the 
hollow  cathode  neutralizer  relative  to  the  body  of 
the  thruster  and  is  often  referred  to  as  the  cathode 
plane  in  the  literature  [17]. 
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The  length  of  the  acceleration  region  for 
each  case  shown  is  invariant  at  30  mm.  Therefore, 
increases  in  the  anode  potential  result  in  linearly 
increased  electric  fields  within  the  thruster  accelera¬ 
tion  channel.  The  initial  acceleration  begins  10  mm 
within  the  thruster  where  the  magnetic  field  has  a 
value  of  approximately  85%  of  the  centerline  (D  = 
0)  maximum.  The  propellant  acceleration  is  com¬ 
pleted  20  mm  beyond  the  exit  plane  when  the  mag¬ 
netic  field  has  a  value  of  approximately  25%  of  the 
centerline  maximum. 

Significant  acceleration  occurs  outside  the 
Hall  thruster.  Taking  into  account  the  uncertainties 
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Fig.  8.  Axial  ionic  velocity  measurements  at  D  = 
0  mm  for  4  discharge  voltage  conditions. 
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inherent  to  the  velocity  measurements  (±500  ml s), 
the  velocity  increment  imparted  into  the  propellant 
outside  the  Hall  thruster  is  essentially  constant  with 
an  average  value  of  5,000  m/s.  Only  for  the  case  of 
a  1 00  V  discharge  voltage,  does  the  majority  of  the 
acceleration  occur  outside  the  thruster.  Higher  dis¬ 
charge  voltages  appear  to  have  a  constant  percent¬ 
age  of  the  acceleration  occurring  externally.  It  is 
also  informative  to  examine  the  energy  deposited 
into  the  propellant.  In  the  case  of  a  100  V  discharge 
voltage,  approximately  90%  of  the  energy  is  depos¬ 
ited  into  the  propellant  beyond  the  exit  plane.  For 
all  cases  above  160  V  discharge  voltage,  the  frac¬ 
tion  of  energy  deposition  beyond  the  exit  plane  is 
nearer  to  60%.  Thus,  the  majority  of  the  energy 
deposition  into  the  Hall  thruster  propellant  stream 
occurs  outside  the  thruster  body.  However  since  the 
thrust  is  equivalent  to  the  momentum  flux,  the 
majority  (65%)  of  the  thrust  is  still  generated  within 
the  thruster  body  in  all  cases  but  that  of  the  100  V 
discharge  voltage.  In  all  the  cases,  approximately 
60  eV  appear  to  not  contribute  to  the  propellant 
acceleration.  This  value  is  constant  to  within  the 
uncertainties  of  the  velocity  measurements  and 
implies  that  the  mechanism  responsible  for  this  loss 
is  invariant  with  the  applied  anode  potential.  This 
energy  loss  is  therefore  most  likely  a  product  of  the 
anode  and  cathode  potential  falls  and  other  mecha¬ 
nisms  inherent  to  the  thruster  design.  It  must  be 
noted  that  the  measurements  presented  thus  far  are 
limited  to  axial  velocities  and  can  not  account  for 
losses  due  to  plume  divergence. 

Several  radial  profiles  of  the  axial  ionic 
velocity  in  Fig.  9  illustrate  the  radial  variation  of  the 
measured  axial  velocities  for  a  discharge  voltage  of 
160  V  at  two  locations  in  the  plume  and  within  the 
thruster  for  a  discharge  voltage  of  200  V.  The  width 
of  the  acceleration  channel  is  approximately  12  mm 
(-6  mm  <  D  <  6  mm)  and  it  is  quite  evident  that  the 
plume  is  sufficiently  diverged  at  the  axial  locations 
examined  such  that  the  axial  velocity  profile  is  flat 
which  strongly  implies  lines  of  constant  potential  in 
the  radial  direction. 

Radial  velocity  measurements  were  per¬ 
formed  in  the  plume.  Combined  with  the  above 
axial  velocity  measurements  with  the  knowledge 
that  the  axial  velocity  appears  to  be  independent  of 
D,  vector  plots  of  the  near  plume  may  be  con¬ 
structed.  Figure  10  shows  a  vector  plot  constructed 
for  a  200  V  discharge  voltage.  In  this  case,  radial 
velocities  vary  linearly  with  D  with  near  zero  veloc¬ 
ity  at  D  =  0  mm  and  peak  at  values  above  6,000  m/s 
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Fig.  9.  Profiles  of  the  axial  ionic  velocity  for  sev¬ 
eral  discharge  voltages  and  locations. 

as  close  as  Z  =  1 3  mm  and  D  =  8  mm.  A  strong  fea¬ 
ture  of  a  Hall  thruster  plume  is  the  central  core  of 
the  plume.  The  propellant  stream  exits  the  thruster 
in  an  annulus,  but  a  intense,  optically  emitting  coni¬ 
cal  feature  starting  near  the  central  magnetic  core 
and  extending  a  distance  into  the  vacuum  chamber 
is  especially  evident  at  higher  discharge  voltages 
and  extends  several  meters  in  vacuum  facilities  with 
lower  background  pressures  than  those  available  at 
Stanford  University.  This  inward  focus  of  the  diver¬ 
gent  propellant  flow  likely  accounts  for  this  com¬ 
monly  observed  plume  structure. 

Ionic  Line  Shape  Analysis 

The  hyperfine  splitting  constants  for  the 
ionic  xenon  5d[4]7/2-6p[3]5/2  transition  at  834.7 
nm  are  only  known  for  the  upper  6p[3]5/2  state. 
Similarly,  the  isotope  shifts  for  this  transition  are 


Fig.  10.  Estimated  flow  vectors  from  radial 
velocity  measurements  for  a  200  V  discharge 


unknown.  In  addition  to  the  lack  of  spectroscopic 
information,  the  plasma  environment  introduces 
uncertainties  in  the  determination  of  an  ionic 
kinetic  temperature  from  an  unsaturatcd  line  shape. 
First,  the  time  averaged  broadening  of  the  excita¬ 
tion  spectra  can  not  be  due  to  plasma  oscillations. 
King  has  shown  in  energy  analyzer  studies  that  the 
axial  velocity  of  the  exiting  ions  have  an  energy  dis¬ 
tribution  of  approximately  10  eV  due  to  the  plasma 
oscillations  within  the  Hall  thruster  [26].  Second, 
the  low  density  of  the  plasma  may  preclude  a  Max¬ 
wellian  velocity  distribution  among  the  ions.  If  so, 
the  concept  of  a  kinetic  temperature  is  not  valid. 
The  issue  of  the  distribution  of  ionic  velocities  is 
minimized  by  examining  the  fluorescence  spectra  in 
the  radial  direction.  The  spectra  is  taken  from  the 
position  with  the  minimum  measured  velocity, 
approximately  100  m/s  at  a  location  of  D  =  0  mm 
and  Z  =  13  mm.  The  ion  population  is  assumed  to 
be  Maxwellian,  or  at  least  frozen  into  a  close  fac¬ 
simile. 

Despite  the  inherent  limitations,  an  esti¬ 
mate  of  the  kinetic  temperature  is  still  possible  if 
uncertainties  of  40-70%  are  acceptable  |27J.  An 
unsaturated  fluorescence  trace  is  shown  in  Fig.  1 1 
and  compared  to  a  model  developed  by  Cedolin 
[28].  The  model  uses  the  5d[3]7/2  lower  level 


Fig.  11.  Model  fit  to  unsaturated  fluorescence 
trace  at  200  V,  D  =  0  mm,  and  Z  =  X  mm. 
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hyperfine  spin  splitting  and  isotopic  shift  data  from 
the  605.1  nm  transition  and  the  measured  splitting 
data  for  the  upper  6p[3]5/2  level.  Lorentzian  broad¬ 
ening  is  neglected  and  only  Doppler  broadening  is 
considered.  The  best  fit  of  this  model  predicts  a 
kinetic  temperature  of  approximately  450  K.  The 
model  still  does  not  completely  predict  the  outlying 
features,  but  this  is  expected  since  the  incorrect  spin 
splitting  constants  are  used  for  the  lower  level.  As  a 
final  note  on  the  uncertainty  of  this  model,  if  hyper¬ 
fine  splitting  is  ignored  and  only  the  isotope  shifts 
corresponding  to  the  values  for  the  605. 1  nm  transi¬ 
tion  are  used,  the  model  predicts  a  kinetic  tempera¬ 
ture  of  approximately  750  K. 

The  kinetic  temperature  of  the  ions  is 
between  450  and  750  K.  The  uncertainty  of  this 
measurement  is  in  large  part  due  to  the  uncertainties 
of  the  spectral  data  as  well  as  due  to  the  noise  in  the 
fluorescence  signal.  A  similar  measurement  in  the 
plume  of  a  SPT-100  by  Manzella  yielded  a  kinetic 
temperature  of  approximately  800  K  (17].  It  should 
be  noted  that  Manzella  used  the  incorrect  value  of  J 
for  the  lower  state  which  was  first  misidentified  by 
Humphreys  |29|  and  propagated  by  Moore  [30] 
before  finally  being  corrected  by  Hansen  and  Pers- 
son  1 18].  This  may  account  for  the  slightly  higher 
kinetic  temperature;  although,  the  uncertainty  is 
such  that  the  difference  in  the  temperatures  is  prob¬ 
ably  inconsequential. 

Neutral  LIF  Velocimetry 

Figure  12  shows  axial  neutral  velocity 
measurements  within  the  acceleration  channel  of 
the  Hall  thruster.  All  four  cases  examined  show 
very  similar  behavior.  The  initial  velocity  near  the 
anode  is  very'  low.  The  neutral  velocity  slowly  rises 
until  a  position  approximately  20  mm  within  the 
thruster.  At  this  point  where  the  ion  acceleration 
also  begins,  the  neutrals  are  accelerated  at  a  higher 
rate  until  near  the  exit  plane  where  the  acceleration 
appears  to  slow  and  even  reverses  when  the  thruster 
is  operated  at  discharge  voltages  of  200  V,  or  above. 
The  decrease  in  the  neutral  velocity  is  best 
explained  by  the  thruster  ingestion  of  background 
neutral  xenon.  Since  the  effect  appears  to  grow  with 
increased  discharge  voltage,  it  is  possible  that  a  por¬ 
tion  of  the  flow  reflected  from  nearby  vacuum  facil¬ 
ity  walls  is  ingested  by  the  thruster.  Once  the 
background  neutral  xenon  is  in  the  vicinity  of  the 
discharge,  the  high  electron  density  near  the  exit 
plane  will  collisionally  excite  the  background 
atoms.  Some  of  these  atoms  will  be  ionized  and  the 
external  electric  fields  produced  by  the  thruster  will 


Z  (mm) 


Fig.  12.  Axial  neutral  velocity  measurements 
at  D  =  0  mm  for  several  conditions. 

accelerate  the  resulting  ions  downstream,  albeit 
with  lower  final  energies  than  ions  created  within 
the  thruster.  Some  of  the  atoms  will  be  electroni¬ 
cally  excited,  but  not  be  ionized  and  will  optically 
decay  to  the  6s[3/2]^  metastable  level  where  they 
could  be  sampled  during  neutral  LIF  velocimetry 
measurements. 

Row  from  the  cathode  can  be  eliminated 
as  the  source  of  the  apparent  slowing  down  of  the 
neutrals  for  several  reasons.  First,  the  flow  from  the 
cathode,  although  15%  of  the  anode  flow,  is  exiting 
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from  a  2  mm  orifice  approximately  12  cm  above  the 
sample  volume.  Flow  from  the  cathode  should  be 
sufficiently  diffused  and  not  affect  the  neutral 
velocity  measurements.  Second,  the  cathode  is 
angled  30°  from  the  front  plate  pointed  down¬ 
stream.  It  is  difficult  to  see  how  flow  from  the  cath¬ 
ode  could  affect  neutral  velocity  measurements 
within  the  acceleration  channel.  It  is  therefore 
almost  certain  that  random  neutral  flux  from  the 
chamber  background  is  responsible  for  the  apparent 
drop  in  neutral  velocity  seen  near  the  exit  plane  in 
Fig.  12. 

Due  to  the  highly  nonequilibrium  nature  of 
the  Hall  thruster,  it  is  important  to  understand  the 
apparent  acceleration  of  the  neutrals.  The  plasma 
within  the  Hall  thruster  is  required  to  be  very  dif¬ 
fuse  by  the  constraint  that  the  magnetic  field 
restrain  the  electron  flux  to  the  anode.  This  require¬ 
ment  appears  to  be  in  force  when  the  neutral  veloc¬ 
ity  measurements  are  compared  with  the  ionic 
velocity  measurements.  The  disparate  velocities  of 
the  ions  and  neutrals  strongly  suggest  that  the  neu¬ 
tral  and  ion  populations  are  not  coupled.  As  such, 
the  apparent  acceleration  of  the  neutrals  may  actu¬ 
ally  be  an  artifact  of  the  time  of  flight  of  the  neutrals 
through  the  volumetric  zone  of  ionization.  Slower 
neutrals,  or  neutrals  that  travel  a  longer  effective 
path  length  due  to  collisions  with  the  walls  of  the 
acceleration  channel,  have  a  proportionally  greater 
chance  being  ionized  than  do  neutrals  in  the  high 
energy  portion  of  the  velocity  distribution.  There¬ 
fore,  neutrals  from  the  high  energy  tail  of  the  veloc¬ 
ity  distribution  are  more  likely  to  reach  the 
upstream  portions  of  the  acceleration  channel.  In 
this  case,  there  is  no  actual  acceleration  of  the  neu¬ 
trals,  but  rather  a  depletion  of  the  slower  moving 
neutrals  by  ionization.  This  has  been  shown  to  be 
possible  by  integration  of  the  Boltzmann  equation 
1311. 

Electric  Field  Calculations 

The  axial  component  of  the  electric  field 
calculated  by  determining  the  ion  kinetic  energy 
from  the  velocimetry  data  in  Fig.  8  and  differentiat¬ 
ing  to  produce  an  effective  electric  field  as  shown  in 
Fig.  13.  It  is  compared  to  electric  field  data  calcu¬ 
lated  from  emissive  probe  plasma  potential  mea¬ 
surements  |8).  The  two  data  sets  show  remarkable 
similarity  indicating  that  the  ionization  and  acceler¬ 
ation  regions  within  the  thruster  are  separated  to 
some  degree.  The  differences  between  the  two  sets 
of  curves  are  due  to  the  differences  between  the 
data.  Ion  creation  occurs  in  a  volume  and  therefore 


the  LIF  derived  plasma  potentials  are  only  valid 
beyond  the  region  of  ion  creation  which  appears  to 
lie  near  Z  =  -10  mm.  Unfortunately,  there  is  proba¬ 
bly  some  ionization  occurring  throughout  the  accel¬ 
eration  channel  which  accounts  for  some  of  the 
differences.  Some  ions  are  likely  created  in  regions 
where  the  electric  field  is  high  and  this  masks  the 
detection  of  local  electric  fields  using  velocimetry 
data,  especially  if  the  zone  of  ion  creation  is  large 
and  the  ion  velocities  are  subsequently  distributed. 

Another  issue  with  determining  plasma 
potentials  and  electric  fields  from  velocimetry  data 
is  that  the  uncertainties  associated  with  the  ionic 
velocity  measurements  are  magnified  by  the  deriva¬ 
tion  of  potential  and  subsequent  electric  field  data. 
With  the  plasma  potential  probe,  the  primary  issue 
is  the  degree  of  disturbance  by  the  probe  while  the 
LIF  velocimetry  is  non-intrusive,  or  less  intrusive 
with  the  addition  of  the  slot  in  the  insulator. 

Both  data  sets  show  similar  trends.  Some 
of  the  features  are  different,  particularly  those  more 
than  10  mm  within  the  Hall  thruster  exit  plane. 
Between  Z  =  -10  mm  and  the  anode,  there  exists 
only  a  small  population  of  ions  with  low  axial 
velocities  with  uncertain  creation  and  loss  rates. 
Velocimetry  studies  prove  inadequate  to  study  the 
potential  and  electric  field  in  this  region,  but  poten¬ 
tial  probe  measurements,  assuming  that  they  do  not 
significantly  perturb  thruster  operation,  provide  bet¬ 
ter  measurement  of  the  plasma  potential  and  subse¬ 
quently  a  better  indication  of  the  electric  field.  It  is 
believed  that  the  probe,  undermeasures  the  higher 
gradients  of  the  plasma  potential  due  to  the  probe 
body  perturbing  the  plasma  flow  field.  This  would 
explain  the  why  the  electric  fields  measured  by  the 
probe  are  somewhat  less  than  those  determined 
from  velocimetry  studies. 

It  is  interesting  to  note  that  the  trends  of 
both  measurements  agree,  particularly  near  the  exit 
plane  where  a  very  low  electric  field  is  indicated. 
This  sudden  change  in  the  electric  field  at  the  exit 
plane  is  visible  in  both  data  sets  and  increases  with 
increasing  discharge  voltage.  If  it  can  be  assumed 
that  the  axial  current  is  a  constant  and  since  the 
local  radial  magnetic  field  is  continuous,  the  dips  in 
the  electric  field  can  be  attributed  to  changes  in  the 
local  plasma  conductivity  indicating  either  a  sharp 
drop  in  the  electron  temperature,  or  a  rise  in  the 
electron  density.  The  measured  electron  tempera¬ 
ture  does  not  support  the  former;  therefore,  a  sharp 
rise  in  the  electron  density  may  be  responsible.  The 
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Fig.  13.  Comparison  of  axial  electric  field  calcu¬ 
lated  from  plasma  potential  and  LIF  velocimetry 
data. 

most  likely  explanation  for  a  rise  in  electron  density 
is  that  due  to  the  high  back  pressure  within  the  vac¬ 
uum  facility  (ICC1  Torr)  there  may  be  a  second  ion¬ 
ization  zone  near  the  exit  plane  of  the  thruster.  In 
flowing  electrons  ionizing  background  neutrals 
could  produce  a  higher  plasma  density  in  this 
region  thus  increasing  the  local  plasma  conductivity 
and  lowering  the  local  electric  field  as  seen  in  Fig. 
13.  Alternatively,  there  may  exist  a  region  near  the 
exit  plane  of  increased  plasma  turbulence  which 
may  produce  increased  plasma  conductivity. 


Conclusions 

Measurements  of  ion  and  neutral  veloci¬ 
ties  were  performed  in  the  plume  and  into  the  inte¬ 
rior  of  the  thruster  through  a  1  mm  wide  slot  in  the 
outer  insulator  wall.  From  these  measurements, 
information  on  propellant  energy  deposition,  elec¬ 
tric  field  strength,  and  flow  divergence  were 
extracted. 

Ionic  velocity  measurements  of  axial 
velocity  both  inside  and  outside  the  thruster  as  well 
as  radial  velocity  measurements  outside  the 
thruster,  were  performed  using  LIF  with  nonreso¬ 
nant  signal  detection  using  the  ionic  xenon 
5d[4]7/2  -  6p[3]s/2  excitation  transition  while 
monitoring  signal  from  the  6s[2]3/2-6p[3]5/2 
transition.  Neutral  velocity  measurements  were 
similarly  performed  in  the  interior  of  the  Hall 
thruster  using  the  6s[3/2]£  -6p[3/2]2  transition 
with  resonance  fluorescence  collection.  Most  veloc¬ 
ity  measurements  used  saturated  transitions  to 
improve  the  signal  to  noise  ratio.  One  radial  trace  of 
the  ionic  transition  was  taken  in  the  linear  fluores¬ 
cence  region  and  yielded  an  ionic  translational  tem¬ 
perature  between  450  and  750  K.  However,  since 
the  hyperfine  structure  constants  are  not  known  for 
the  5d[4]7/2  level,  the  constants  for  the  5d[3]7/2 
level  were  used  instead.  This  result  should  therefore 
be  viewed  with  caution.  The  ionic  velocity  data 
allowed  for  the  measurement  of  the  energy  depos¬ 
ited  into  the  propellant  stream  which  was  calculated 
as  the  kinetic  energy  of  the  ions.  An  effective  elec¬ 
tric  field  was  calculated  from  the  ionic  kinetic 
energy.  Electric  field  values  derived  from  the  ion 
velocity  measurements  show  that  the  fields  are  rela¬ 
tively  smooth  and  peak  near  the  exit  plane.  Maxi¬ 
mum  field  strengths  were  found  to  be  approach 
8,000  V/m.  At  the  exit  plane,  the  axial  electric  fields 
show  a  dip  for  nearly  all  the  discharge  voltages 
examined. 
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